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Abstract 

We investigate the back reaction of cosmological perturbations on the evolution of the Universe 
using the renormalization group method. Starting from the second order perturbed Einstein's 
equation, we renormalize a scale factor of the Universe and derive the evolution equation for the 
effective scale factor which includes back reaction due to inhomogeneities of the Universe. The 
resulting equation has the same form as the standard Friedman-Robertson- Walker equation with 
the effective energy density and pressure which represent the back reaction effect. 
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I. INTRODUCTION 



Owing to the nonlinear nature of Einstein's equation, fluctuations of the metric affect the 
evolution of the background space time. In cosmology, this effect is expected to be important 
when we consider the evolution of large scale nonlinear structures and the evolution of the 
early Universe. This is cosmological the back reaction problem and has been studied by 
several authors |, |, |, |, §, §, 0, |, |, |10|, |TT|. 



One of the main difficulty in the cosmological back reaction problem is connected with 
the gauge freedom of perturbation variables. The back reaction effect appears from the 
second order quantities in a perturbation expansion and we must use the second order gauge 
invariant quantities to avoid the gauge ambiguity of the back reaction problem. Abramo 
and co-workers || |7|, [| derived the gauge-invariant effective energy momentum tensor of 
cosmological perturbations and applied their formalism to the inflationary universe. They 
discussed the effect of the inhomogeneity on the background Friedman-Robertson- Walker 
(FRW) universe, but did not derive solutions of an effective scale factor for the FRW universe 
with the back reaction. 

In our previous papersJlO], [□]], the renormalization group method [[L2|, [3~3| , jlij |i~5H was 



applied to the cosmological back reaction problem. We start from the following perturbation 
expansion of the metric 

(0) (1) (2) 

9ab = 9 ab + 9 ab + 9 ab + i I 1 ) 

(0) 

where g ab is the background FRW metric and represents a homogeneous and isotropic space. 
V ab is the metric of the first order linear perturbation with a fc) = where (• • ■) means the 

(2) 

spatial average with respect to the background FRW metric, g ab is the second order metric 
and this part contains nonlinear effects caused by the first order linear perturbation. This 
nonlinearity produces homogeneous and isotropic zero modes as part of the second order 
metric. That is, 

(fab) + 0. (2) 

As the zero mode part of the metric has the same symmetry as the background FRW 
universe, it must be interpreted as a part of FRW metric. Hence we redefine the background 
metric as follows: 

(B) (0) (2) 

9 ab = 9ab + {dab/- ( 3 ) 

(B) 

We used the renormalization group method to define the new background metric g ab which 



includes the back reaction effect. In previous papers |TU| |TTf, we first solved perturbation 



equations and obtained g ab and g ab . Then the renormalized effective scale factor with the 
back reaction effect was obtained. 

But it is more convenient to obtain the renormalized Einstein's equation for the zero mode 
variables from the first to investigate the back reaction effect. This is equivalent to average 
Einstein's equation and derive the evolution equation for the effective FRW universe with 
the back reaction effect. In this paper, we apply the renormalization group method directly 
to the perturbed Einstein's equation. We does not use an explicit form of the first and the 
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second order solutions of perturbations. All that we need is how constants of integration 
enter in the solution of perturbations. 

The plan of this paper is as follows. In Sec. II, we introduce the gauge ready formalism 
of cosmological perturbations and derive the second order Einstein's equation for gauge 
invariant variables. In Sec. Ill, the renormalization group method is applied to the zero 
mode part of the Einstein's equation and the evolution equation for the renormalized scale 
factor is derived. In Sec. IV, we investigate the solution of the back reaction equation. 
Sec. V is devoted for summary and discussion. We use units in which c = h = 8nG = 1 
throughout the paper. 



II. THE SECOND ORDER PERTURBATION FOR THE UNIVERSE WITH 
IDEAL FLUID 



To circumvent the gauge ambiguity of cosmological perturbations and to obtain a 
gauge independent interpretation of the cosmological back reaction problem, we must 
use gauge invariant description of cosmological perturbations [16, [17]]. Abramo and co- 
workers obtained the gauge invariant effective energy-momentum tensor for cosmological 
perturbations || |7|, |[ [| which is quadratic with respect to the first order variables. To 
apply the renormalization group method to the back reaction problem, we must obtain the 
second order zero mode perturbation in a gauge invariant manner [ITTJL [TT||. In this section, 



we introduce the gauge ready formalism ||T8[ and derive the evolution equation for the second 
order zero mode perturbation which is invariant under the first order gauge transformation. 



A. the gauge ready form of cosmological perturbations 

We consider a spatially flat FRW universe with perfect fluid of which equation of state 
is given by p = w p, where w is assumed to be constant. The background scale factor of the 
Universe and the energy density are given by 

«W = cMa, 9=-^, (4) 

where C and Cq are constants of integration. The energy density of the fluid p is obtained 
from the zeroth order of the following mass conservation equation 

piT^ v /37J M a j = o, (5) 



where u a is four velocity of the fluid. We restrict our attention only to the scalar type 
perturbation and the metric of the first order perturbation can be written as 

(i) (-2<f> a(t)B, \ 

a {t)B ti a 2 (t) [(1 - 2^) 5 i3 + 2E >ij ) J ' W 

Let us consider the infinitesimal coordinate transformation 
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x m = x o + £0 x n = ^ + 6 t 3 £ / 7 ) 

a z 
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Under this transformation, the perturbation variables receive the following gauge transfor- 
mations: 



0' = 0-£°, 

f = v + m°, 



B' — B — a 



+ 



E' = E 



The comoving three velocity u$ = v » of the fluid transforms as 

f ' = t> + a 2 



.2f e^' 



and the velocity potential defined by x = t> + aB transforms as 



(8) 



(9) 



(10) 



From these transformations, we can make the following gauge invariant combination of 
perturbation variables: 



$ = <p + a(B- aE) 
X = x-a(B - aE), 
Now the vector 



V = ip-aH(B-aE), 



e c = (-x, s*E d ) 



(11) 



(12) 



transforms as four vector under the coordinate transformation. The gauge transformation 
induced by the coordinate transformation using the vector (O) is 



(p' = (p + X = ® + X, 
E' = E - E = 0, 

x' = x - x = o. 



a a 



(13) 



Hence in this gauge, each components of the metric correspond to the gauge invariant 
variables and a comoving gauge condition x = is satisfied. The choice of the vector ( |T2"D 
is not unique. The vector 



tt= {-a(B-aE), P'Ej) (14) 
also transforms as a four vector and the gauge transformation using this vector is 



+ 



a(B - aE) 



ijj' = ij) - aH(B - aE) = 



B' = 0, 

x' = x-a(B- aE) = X. 



E' = 0, 



(15) 
(16) 



This gauge corresponds to the longitudinal gauge and components of the metric coincide 
with gauge invariant variables $ and This gauge is often used to investigate cosmo logical 
perturbations [17[], and Abramo and co-worker || [7|, ||, [| also used this gauge to evaluate 
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the gauge invariant energy momentum tensor of cosmological perturbations. But for the 
purpose of the back reaction problem, this gauge is not suitable because the contribution of 
long wavelength limit mode of the perturbations to the effective energy momentum tensor 
does not vanish automatically . The perturbation of which wavelength is infinite cannot 
be distinguished from the background and has nothing to do with the back reaction effect. 
But in the longitudinal gauge, we have contribution of the long wavelength limit perturbation 
to the second order quantities. We must subtract this to obtain correct back reaction effect. 
On the other hand, in comoving gauge, the second order perturbation becomes zero in the 
long wavelength limit. In this paper, therefore, we adopt the gauge ready form in the 
comoving gauge to investigate the back reaction problem. 

To simplify the calculation, we choose B = and the metric up to the second order in 
the gauge (|13~D is given by 



ds 2 = -(1 + 20 + 20 2 ) dt 2 - 2a 2 E ;i dx l dt + a 2 (t)(l -2ip- 2ip 2 ) dx 2 . (17) 

As the each component of the first order metric perturbation corresponds to the gauge 
invariant variable in this gauge, the second order zero mode variables 02 and ip 2 determined 
by the first order quantities are also invariant under the first order gauge transformation up 
to the second order. 

Using the metric flT7j ) with the comoving condition \ = 0, the first order Einstein equa- 
tions are 

1 V 2 ^ + HV 2 E = \ (18a) 
ij + H(p = 0, (18b) 
$ + 3Hip + H<j>- 3wH 2 (f) = ^ p, (18c) 

E + 3HE + \(^-(f)) = 0. (18d) 
a 2 

We can obtain the following evolution equation for the spatial curvature perturbation ip 

ip + 3H^j - — V 2 V> = 0. (19) 
a 2 

B. the evolution equations for the second order zero mode perturbation 

The second order Einstein equations for the zero mode variables 02 and ip2 are 

6H(^ 2 + H<p 2 ) = -(G° ( ( g)) - p, (20a) 

1 (2) 

2[^ + #(3</>2 + 2 ) - 3H 2 w<p 2 ] = -~(G\(g)) + w p, (20b) 



where 

(2) m / ^ 
(G° (V)) = ( -3^ 2 + 12H^ - ^V 2 ij - 6HijV 2 E 



a 2 



- 3H 2 a 2 EV 2 E - 2F0 (- V 2 E + Qtjjj - 12F 2 2 ) , (21) 
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(G'jCg)) = ( -^V 2 ^ + ^ + 4^# + 3ff0) - 2V 2 E (3H1P + 

2^V 2 E - 2a 2 HEV 2 E - (2H + 5H 2 )a 2 EV 2 E + -£V 2 - 20'0 

3 



+ 3HiP - ^V 2 ^ - i (V 2 £ + 2#V 2 £) + 2H<pj 



^0V 2 - 40 2 (3^ 2 + 2H )\ 5V (22) 



3a 2 

Eq.( [20a| ) is a time-time component and Eq. (|20b|) is a space-space component of the spatially 
averaged Einstein's equation (G a b) = (T a b)- The spatial average is defined by 

1 r„ 



(A) = — I d 6 xA } 



where V is the volume of a sufficiently large compact domain and periodic boundary condi- 
tions for perturbations are assumed. Tl 
the mass conservation (pf) and given by 



tions for perturbations are assumed. The second order energy density p is determined by 



(2) ft I \ {0) (?l i \ I OZ/WV , W ( ( P 2 ) , /l , /nn\ 

/9 = (1 + w) p (3^ 2 - 2 ) + 3(V> p ) + — — -— — + ca(l + w) p , (23) 

v ; p 



where C2 is a constant. 



III. THE RENORMALIZATION AND THE EFFECTIVE EINSTEIN'S EQUA- 
TION 

In this section, we derive the effective Einstein's equation which describes the evolution 
of the effective FRW universe with the back reaction effect using the renormalization group 
method. The spatially averaged Einstein's equation up to the second order is 



(0) 



- 3H 2 + p + QH(ip 2 + Hfa) + (1 + w) p (3^ 2 - 02 + c 2 ) = -Pbr 



(0) 



- 2H - 3H 2 - w p + 2 ^2 + H(3ip 2 + <h) ~ 3H 2 w4> 2 

- w(l + w)p(3ip 2 - 02 + c 2 ) = Per, 



(24a) 



(24b) 



where 



( 2 )n (1) (1) 



w 



/ (1) 2\ 
(P ) 



2(1 + w) (o) 



2 ( ( P )2 ) 



v 1 p 



(25a) 
(25b) 
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By introducing a new time variable dr = (1 + 4> 2 ) dt, the Einstein equations (|24a ) and fl24b|) 
become 



3H 2 + 6Hijj 2 + ( p [1 + (1 + w) (3^2 - h + c 2 )] = -Per, (26a) 
2ij - 3# 2 + 2 (^ 2 + 3H^ 2 ) - w ( p [1 + (1 + w) (3^ - 02 + c 2 )] = Per, (26b) 



where ' = 4- . 

dr 

The spatially averaged line element up to the second order perturbation is 

(ds 2 ) = -c/r 2 + a 2 (r)[l-2^ 2 (r) + 2^ 2 (r )] dx 2 , a{r) = Cr^i) , (27) 
where ro is a time at which the initial condition is set, and it is always possible to write 

(2) 

the zero mode metric as Eq.(|2T|) by choosing a constant of integration c 2 contained in p as 
C2 = — 3ip 2 (T ). From this form of the line element, we can observe that the effective scale 
factor for the FRW universe is 

flcff(r) =rmc [1-Mt) + Mto)]- ( 28 ) 

We regard the second order perturbation ip 2 as a secular term and apply the renormalization 
group method [|T2"|, |T^, [T4|, IISJ to absorb it to the zeroth order constant of integration C. We 
redefine the zeroth order integration constant C as 

C = C(ji) + 6C(n;r ), 

where p is a renormalization point and 5C is a counterterm which absorbs the secular 
divergence of the solution. We assume that SC is the second order quantity The effective 
scale factor (|28|) up to the second order of the perturbation can be written 

aeff(r) = (C(p) + 8C(p; r ))(l - ^ 2 {t) + ^ 2 {p) ~ Mv) + M T o)) 

= C(p)(l - i) 2 {r) + fo{n)), 

where we have chosen the counterterm 5C so as to absorb the (ip 2 (/■*) — ^(fo^-dependent 
term: 

8C(p; t ) + C{p){-^ 2 {p) + ^(r )) = 0. (29) 

This defines the renormalization transformation 

7V ro : C(r ) 1— C(ji) = C(r ) - C(/i) - ^(r )] , (30) 

and this transformation forms Lie group. The renormalization group equation is obtained 
by differentiating the both side of the transformation with respect to p and setting r = p: 

dC^ = -C(p)^, (31) 
dp dp 

and its solution is 

C(p) = C{t ) e -^M+<Mro)_ (32) 
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By equating p, = r in the effective scale factor, the renormalized scale factor is given by 

a R ( T ) = r^T C(t) = a(r) e -^(r)+^(r )_ (33) 



By substituting the relation a(r) = aR(r)e^ T ' ^ 2<T °) into ( 26a|) and (|26b| ) , and keeping 



terms up to the second order of the perturbation, we obtain the following equations for the 
renormalized scale factor or: 

- 3H 2 + p = -Per, (34a) 
-2H- 3H 2 - w Po = per, (34b) 

where H = ^ and 

PO = [1 - (1 + v)M jj^y = [(1 + ^ a?i]1+ .- 

( P4a| ) and ( f34b| ) are the main result of this paper. The second order curvature variable ^2 
disappears in the evolution equation by the procedure of the renormalization. The second 
order lapse function 02 remains in the expression po, but this variable corresponds to the 
second order gauge freedom to parameterize time and we can freely choose the form of this 
function. Eqs.( |34a| ) and ( |34b| ) have the same form as the standard FRW equations. On the 
right-hand side, the back reaction effect appears as source terms pbr and pbr, and their 
explicit form is determined by solving the first order perturbation. 



IV. SOLUTIONS OF THE EFFECTIVE EINSTEIN'S EQUATION 



In this section, we solve the effective Einstein equations ( |34a| ) and ( |34b|) , and investigate 



how the inhomogeneity affects the expansion of the Universe. By eliminating the variable 
Po, we have 

-2H -3(1 + w)H 2 = wp BR + p BR . (35) 

For w 7^ — 1, by using a new variable x = (a^) U 1+w \ 

3 

x = --(1 + w)(wp BR + p BR )x. (36) 
We solve this equation for various value of w. 

A. vacuum energy case w = — 1 

The universe expands exponentially a oc e Ht ,H = const., and we must treat this case 
separately. In this case, the perturbation of energy density is zero and the effective Einstein's 
equation becomes 

-3H 2 + p = -p BR = -(G°o( {1 g ) )), 

-2H- 3H 2 + Po = p BR = -£<G\(0)>, Po = const. (37) 
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For long wavelength perturbation k <C aH, the first order growing mode solution is given 

by 

k k 

where C k is a constant in time. The source terms of the effective Einstein's equation are 

Pbr ~ \ £ ^l^l 2 ' Pbr ~ £ ^l^l 2 = - ^Pbr- 

This is equivalent to a negative spatial curvature term. For the slow roll inflationary phase 
driven by a scalar field, it was shown that the back reaction effect by the long wavelength 
perturbation is equivalent to positive spatial curvature [PH . For vacuum energy case, fluc- 
tuation of the energy density becomes identically zero and this leads to the different back 
reaction effect compared to the inflation with a scalar field. 

For short wavelength mode k 3> aH, the first order curvature perturbation is given by 



w, ~ -Ji cos h ( -A- ) e lk x 
^fka \aH) 



and 

7 v— \ Cfe 7 v-^ |Cfe| 

k R k R 

The back reaction effect is equivalent to a radiation fluid and this result is the same as the 
inflation with a scalar field. 



B. dust case w = 

The exact solution of the first order growing mode is 

^Vae' h = 0, E = - 3 -r-^YC k e* k *, $ = -f £ % C k e ik *, (38) 
fc fc fc 

where C k is a constant in time. Using this solution, the source terms of the effective Einstein 
equation are 

13^fc 2 |C fc | 2 13 ^k 2 \C k \ 2 1 

^--25^"^"' PBR -3^L^-l m ' (39) 
fc fc 

This is the same equation of state as a positive spatial curvature term. The equation ( |36|) 
becomes 

x = -— x~ l '\ A=—Yk 2 \C k \\ (40) 

4 ' 3-25^ 1 1 K ' 

k 

and the solution is 

IE. . AE 

ttR = 9A^ -cosr]), t = g^iv-smvh ( 41 ) 

where E is a constant of integration. This solution is the same as a closed FRW universe 
with dust. Therefore, the inhomogeneity works as a positive spatial curvature in the dust 



dominated universeH, 10| and the Universe will recollapse. 
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C. w/0 case 

For long wavelength perturbations k <C aH, using the long wavelength expansion, the 
first order solution up to 0(k 2 ) is given by 

^ 3w + 5 



k 



p ~-l^T^^ Cfce ' (42) 



and 



Pbr 



9w 2 + 30w + 13 ^ A; 2 . 

> — ^ <^fc 



( 3u , + 5)2 ^ a R 2 
162w 3 + 423w 2 + 276w - 13 v-^ k 2 ,„ . , 



3(3w + 5) 2 fe a ^ 2 



The effective Einstein equation (RO) becomes 



3 3m-l 

£ = -(! + w)(wA — 5) x 3 (™ +1 ) , (44) 



where 



yt = _ 9u,» + 3( to +13 E B = _ 162^ + W + 240u,-13 E 

(3u> + 5)2 z —' 3(3u> + 5)2 ^ 

By integrating the equation with respect to r, we obtain 

x 2 . . ^ 3(w + l) 2 (135w 3 + 333«; 2 + 201w;-13)^ ;2l ^ l2 W) 

2 cttV ; ' 8 (3w + l)(3w + 5) 2 * 

(46) 

where i£ is a constant of integration. If there is no inhomogeneity, V e g = and the solution of 
the effective Einstein's equation reduces to that of the background FRW solution. The con- 
dition E > is necessary to reproduce the background the solution when the inhomogeneity 
vanishes. We can obtain qualitative behavior of the effective scale factor by observing the 
shape of the effective potential V e g(x) (FIG. 1). We can summarize the back reaction effect 
of the long wavelength perturbation on the expansion of the Universe as follows: 

• — 1 < w < — 1/3: the expansion of the Universe is accelerating. The inhomogeneity 
reduces the expansion rate. 

• w = — 1/3 : V e g <x In x. The expansion of the Universe is decelerating. The inhomo- 
geneity reduces the expansion rate and the universe will recollapse due to the back 
reaction effect. 

• — 1/3 < w < w* : io* is a real root of the equation 135w 3 + 333u; 2 + 201w — 13 = and 
w* w 0.06. The expansion of the Universe is decelerating. The inhomogeneity reduces 
the expansion rate and the universe will recollapse. 
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FIG. 1: Effective potentials for the long wavelength perturbations, (a) is for — 1 < w < —1/3, (b) 
is for — 1/3 < w < w* and (c) is for w* < w < 1. For — 1 < w < w*, inhomogeneity reduces the 
expansion rate of the Universe. For w* < w < 1, inhomogeneity enhances the expansion of the 
Universe. For w = w*, there is no back reaction effect. 

• w = w* : Vpff = and there is no back reaction. 



• w* < w < 1 : the expansion of the Universe is decelerating. The inhomogeneity 
increases the expansion rate. 

For short wavelength perturbations k 3> a H, WKB type solution for the curvature 
perturbation if> is given by 



1> 



Ch cos(av^ / — ) < ik - 



k\fw a 



a J 



for w > 0, 



Y cosh ( [ —) e ikx for w < 0, 

ky/-w a V J a J 



(47) 



and using this solution, we have 

. v 5-9iu 



Per ~ sign (w ) 



E 



\n 12 



n 12 



, / — , PBR~signw— > 

fc fc 



The equation of state for the back reaction term becomes 



(48) 



(49) 



which is independent of the value w. For w < 5/9, per is positive and this is equivalent to 
a radiation fluid. The effective Einstein's equation becomes 



x 2 



— + V cS (x) = E, V eS 



-sign(w, 



3 (w + 1) 2 (5 - 9w) 



x 2(3w-l) 

2j|C fc | 2 o;WTr, £>0. (50) 
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The back reaction effect of the short wavelength perturbation on the expansion of the 
Universe is summarized as follows (FIG. 2): 
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FIG. 2: Effective potentials for the short wavelength perturbations, (a) is for — 1 < w < 1/3, (b) 
is for 1/3 < w < 5/9 and (c) is for 5/9 < w < 1. For — 1 < w < 1/3,5/9 < w < 1, inhomogeneity 
reduces the expansion rate of the Universe. For 1/3 < w < 5/9, inhomogeneity enhances the 
expansion of the Universe. For w = 1/3, 5/9, there is no back reaction effect. 

• — 1 < w < 1/3 : The inhomogeneity reduces the expansion rate. 

• w — 1/3 : V c fi = 0. The equation of state of the back reaction term is the same as the 
back ground matter field and we have no back reaction from inhomogeneity. 

• l/3<u><5/9: The inhomogeneity increases the expansion rate. 

• w = 5/9 : V e g = 0. There is no back reaction. 

• 5/9<u><1: The inhomogeneity reduces the expansion rate and the Universe will 
recollapse. 

V. SUMMARY AND DISCUSSION 



In this paper, we applied the renormalization group method to the perturbed Einstein's 
equation and derived the evolution equation for the renormalized scale factor. The renor- 
malized scale factor includes the back reaction effect caused by the inhomogeneity of the first 
order perturbation. The resulting equation has the same form as the standard FRW equa- 
tion with the back reaction effect appears as the effective energy density and pressure on the 
right hand side of the equation. For the long wavelength perturbation, the equation of the 
state of the back reaction terms depends on the value w. We have pbr = — Pbr/3, Pbr > 
for w = —1 and this is equivalent to a negative spatial curvature. For w = 0, we have 
Pbr = — Pbr/3, Pbr < and this is equivalent to a positive spatial curvature. On the other 
hand, for the short wavelength perturbation, the equation of state of the back reaction term 
becomes independent of w and pbr ~ Pbr/3. For w < 5/9, pbr > and the equation of 
state is the same as a radiation fluid. 

We derived the effective FRW equation fl34a|) and ( |34b|) without using the explicit form of 
the background solution and the perturbation solution. Usually, the renormalization group 
method is utilized to obtain the globally valid approximated solution of the differential 
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equation and we must prepare naive perturbative solution of the original equation before 
applying the procedure of renormalization. In this paper, we could obtain the evolution 
equation for the renormalized scale factor without solving the perturbation equations. Once 
the effective FRW equation is obtained, it is possible to investigate the evolution of the 
effective scale factor and the back reaction effect by numerical method. Hence this approach 
is useful for the cosmological model with the scalar field in which case obtaining the solution 
of the perturbation is not easy in general. 
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